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Abstract 

Given positive numbers pi < P2 < ■ ■ ■ < Pn, and a real number r let 
L r be the n x n matrix with its i, j entry equal to (pi — Pj)/(pi — Pj). A 
well-known theorem of C. Loewner says that L r is positive definite when 
0 < r < 1. In contrast, R. Bhatia and J. Holbrook, (Indiana Univ. Math. 

J, 49 (2000) 1153-1173) showed that when 1 < r < 2, the matrix L r has 
only one positive eigenvalue, and made a conjecture about the signatures 
of eigenvalues of L r for other r. That conjecture is proved in this paper. 
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1 Introduction 


Let / be a real-valued C 1 function on (0, oo). Let pi < P 2 < ■ • • < p n be any n 
points in (0,oo). The nx n matrix 


Lf(pi,...,p n ) 


7(ft) ~ fiPjY 


Pi ~ Pj 


i,j =1 


(1) 


is called a Loewner matrix associated with /. It is understood that when i = j , 
the quotient in 0 represents the limiting value f'(pi). Of particular interest to 
us are the functions f(t) =t r ,r£ M. In this case we write L r for Lf(p 1 ,... ,p n ), 
where the roles of n and p \,.. .p n can be inferred from the context. Thus L r is 
the n x n matrix 


Lft -Pj 



Loewner matrices are important in several contexts, of which we mention 
two that led to the present study. (The reader may see Section 4.1 of [12] for an 
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excellent discussion of both these aspects of Loewner matrices.) The function 
/ on (0, oo) induces, via the usual functional calculus, a matrix function f{A) 
on the space of positive definite matrices. Let Df(A) be the Frechet derivative 
of this function.This is a linear map on the space of Hermitian matrices. The 
Daleckii-Krein formula describes the action of this map in terms of Loewner 
matrices. Choose an orthonormal basis in which A = diag(pi,... ,p n ). Then 
the formula says that for every Hermitian A' 

Df(A)(X) = L f (p 1 ,...,p n ) oX, (3) 

where A o B stands for the entrywise product [ aijbij] of A and B. 

The function / is said to be operator monotone on (0,oo) if A > B > 0 
implies f{A ) > f(B). (As usual A > 0 means A is positive semidefinite.) A 
fundamental theorem due to Charles Loewner says that / is operator monotone 
if and only if all Loewner matrices associated with / (for every n and for every 
choice pi ,... ,p n ) are positive semidefinite. Another basic fact, again proved first 
by Loewner, says that f(t) = t r is operator monotone if and only if 0 < r < 1. 
See pQ Chapter V. 

Combining these various facts with some well-kown theorems on positive 
linear maps [2j one can see that if / is operator monotone, then the norm of 
Df(A) obeys the relations 

\\Df(A)\\ = \\Df(A)(I)\\ = \\f'(A)l (4) 

and is therefore readily computable. In particular, for the function f(t) = t r if 
we write DA r for Df(A), then ([4]) gives 

||DA , '|| = ||rA r_1 ||, for 0 < r < 1. (5) 

This was first noted in [3j, and used to derive perturbation bounds for the 
operator absolute value. Then in [3] Bhatia and Sinha showed that the relation 
<[5j) holds also for —oo < r < 0 and for 2 < r < oo but, mysteriously, not for 
1 < r < \/2. The case y/2 < r < 2, left open in this paper, was resolved in [J] 
by Bhatia and Holbrook, who showed that here again the relation ([5| is valid. 

One ingredient of the proof in [4j is their Proposition 2.1 which says that 
when 1 < r < 2, the n x n matrix L r has just one positive eigenvalue. We have 
remarked earlier that when 0 < r < 1, the matrix L r is positive semidefinite 
and therefore, none of its eigenvalues is negative. This contrast as r moves 
from (0,1) to (1,2) is intriguing, and raises the natural question about the 
behaviour of eigenvalues of L r for other values of r. Bhatia and Holbrook )4] 
made a conjecture about this behaviour and established a small part of it: they 
settled the cases r = 1,2,..., n — 1 apart from 0 < r < 1 and 1 < r < 2 already 
mentioned. The main goal of this paper is to prove this conjecture in full. This 
is our Theorem 1.1. 

Let A be an n x n Hermitian matrix. The inertia of A is the triple 
In (A) = (7 t(A),C(A),j/(A)), 
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where it {A) is the number of positive eigenvalues of A, £(A) is the number 
of zero eigenvalues of A, and v(A) the number of negative eigenvalues of A. 
Theorem |1.1| describes the inertia of L r as r varies over R. It is easy to see that 
the inertia of L_ r is the opposite of the inertia of i.e. tt(L_ r ) = v(L r ) and 
v{L_ r ) = tt(L r ). So we confine ourselves to the case r > 0. 

Theorem 1.1. Let pi < P 2 < • • • < p n and r be any positive real numbers and 
let L r be the matrix defined in Q). Then 

(i) L r is singular if and only if r = 1,2,..., n — 1. 

(ii) At the points r = 1,2,... ,n, the inertia of L r is given as follows: 

r = 2k => In (L r ) = (k,n — r,k), 

and 

r = 2k — 1 => In ( L r ) = ( k, n — r, k — 1). 

(iii) If 0 < r < n and r is not an integer, then 

[rj = 2k =>■ In (L r ) = (n — k, 0, k ) 

and 

[_rj = 2k — 1 => In (L r ) = ( k, 0, n — k). 

(iv) If r > n — 1, then In (L r ) = In ( L n ). 

(v) Every nonzero eigenvalue of L r is simple. 

It is helpful to illustrate the theorem by a picture. Figure 1 is a diagram 
of the (scaled) eigenvalues of a 6 x 6 matrix L r when pi are fixed and r varies. 
Some of the eigenvalues are very close to zero. To be able to distinguish between 
them the vertical scale has been expanded. 


We have already mentioned that for 0 < r < 1, statement (iii) of Theorem 
0 follows from Loewner’s theorem, and for 1 < r < 2 it was established in g]. 
The case 2 < r < 3 was accomplished by Bhatia and Sano in 7j. We briefly 
explain this work. 

Let Hi be the space 

H\ = ja: = [x\,..., x n ) : ^ Xj = 0 j . (6) 

An n x n Hermitian matrix A is said to be conditionally positive definite if 
{x,Ax) > 0 for all x £ Hi, and if —A has this property, then we say that A is 
conditionally negative definite. Since dinrHi = n— 1, a nonsingular conditionally 
positive definite matrix which is not positive definite has inertia (n — 1, 0,1). 
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Eigenvalues of L r ; n = 6, 0 < r < 10 



Figure 1: 


In J7i| it was shown that when 1 < r < 2, the matrix L r is nonsingular 
and conditionally negative definite. It follows that In ( L r ) = (1,0, n — 1), a 
fact established earlier in [2j. It was also shown in [7] that when 2 < r < 3, 
the matrix L r is nonsingular and conditionally positive definite. From this it 
follows that In (L r ) = (n — 1,0,1). 

More generally, Bhatia and Sano [7] showed that / on (0, oo) is operator 
convex if and only if all Loewner matrices Lf are conditionally negative defi¬ 
nite. This is a characterisation analogous to Loewner’s for operator monotone 
functions. It is well-known that f{t) = t r is operator convex for 1 < r < 2. 

The proof of Theorem o is given in Section 2. We also indicate how the 
proofs for the parts already given in [4] and (7] can be considerably simplified. 
The inertia of the matrix [(p, + PjY] has been studied by Bhatia and Jain in 
[5]- Some ideas in our proofs are similar to the ones used there. 

2 Proofs and Remarks 

Let X be an n x n nonsingular matrix. The transformation A i—»• X* AX on 
Hermitian matrices is called a congruence. The Sylvester Law of Inertia says 
that 

In (. X*AX ) = In A for all X G GL{n). (7) 
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Let D be the diagonal matrix 


D = diag (pi,... ,p n ). 


Then for every r 


L_ r = -D~ r L r D~ r . 


Hence by Sylvester’s Law 


In L r = (* i ,*2,*3) •<=> In = (*3>*2,* i )- 


( 8 ) 

(9) 

( 10 ) 


Thus all statements about In L r for r > 0 give information about In L_ r as 
well. 

Make the substitution pi = e 2xi , Xi £ M. A simple calculation shows that 


e rXi sinh r(xi — Xj) e rXj 
e Xi sinh(:r, — Xj) e x i 

In other words, 

L r = AL r A, 

where A = diag (e^ -1 ^ 1 ,..., gl’’ -1 )®"), and 

~ sinh r(xi — Xj) 
sinh(a;i — Xj) 


( 11 ) 


( 12 ) 


By Sylvester’s Law In L r = In L r . Several properties of L r can be studied via 
L r , and vice versa. This has been a very effective tool in deriving operator 
inequalities; see, the work of Bhatia and Parthasarathy [B] and that of Hiai and 

Kosaki pnnunma. 

When n = 2 we have 




sinh r(x\ — X2) 
^ sinh(xi — X2) 

sinh r(x 1—X2) 
sinh(aii —X2) 


So detL r = r 2 — ■ Thus detL r is positive for 0 < r < 1, zero for 

r = 1, and negative for r > 1. One eigenvalue of L r is always positive, and 
this shows that the second eigenvalue is positive, zero, or negative depending 
on whether 0<r<l,r = l, or r > 1, respectively. This establishes Theorem 
0 in the simplest case n = 2. 

An interesting corollary can be deduced at this stage. According to the two 
theorems of Loewner mentioned in Section 1, / is operator monotone if and only 
if all Loewner matrices Lf are positive semidefinite, and f(t) = t r is operator 
monotone if and only if 0 < r < 1. Consequently, if r > 1, then there exists an 
n, and positive numbers p\,... ,p n such that the associated Loewner matrix Q 
is not positive definite. We can assert more: 


Proposition 2 . 1 . Let r > 1. Then for every n > 2, and for every choice of 
Pi ,..., p n , the matrix L r defined in © has at least one negative eigenvalue. 
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Proof Consider the 2x2 top left submatrix of L r . This is a Loewner matrix. 
By Theorem o it has one negative eigenvalue. So, by Cauchy’s interlacing 
principle, the n x n matrix L r has at least one negative eigenvalue. □ 

The Sylvester Law has a generalisation that is useful for us. Let n > r, and 
let A be an r x r Hermitian matrix and X an r x n matrix of rank r. Then 


In X*AX = In A + (0,n — r, 0). (13) 

A proof of this may be found in [5]. This permits a simple transparent proof of 
Part (ii) of Theorem |1.1| (This part has already been proved in [4].) When r is 
a positive integer we have 


Lr = [pr 1 +pT 2 Pj + ■ ■ ■+P r j ~ 1 } = W*vw, 

where W is the r x n Vandermonde matrix 


W = 


1 

Pi 


Pi 1 


1 

P2 


P2 1 


1 

Pn 


and V is the r x r antidiagonal matrix with all entries 1 on its sinister diagonal 
and all its other entries equal to 0. If r = 2 k, the matrix V has k of its eigenvalues 
equal to 1, and the other k equal to —1. If r = 2k — 1, then k of its eigenvalues 
are equal to 1, and k — 1 are equal to — 1. So, statement (ii) of Theorem o 
follows from the generalised Sylvester’s Law (13). Next we prove statement (i). 

Let Ci,C 2 ,... ,c„, be real numbers, not all of which are zero. Let / be the 
function on (0, oo) defined as 


ft \ Y^ x ^i 
/O) = l^ c i - - 


i=i 


x - Pj 


(14) 


Theorem 2.2. Let r be a positive real number not equal to 1, 2,..., n— 1. Then 


the function f defined in (If) has at most n — 1 zeros in (0, oo). 


Proof Let rq < ?"2 < • • • < r m , and let ai,..., a m be real numbers not all of 
which are zero. Then the function 

m 

9i x ) = ^2 a jX r \ (15) 

i=i 


has at most m — 1 zeros in (0, oo). This is a well-known fact, and can be found 
in e.g., [15], p.46. 


Now let / be the function defined in (141 and let 


g(x) = f{x) Y[(x- Pj ). 
i=i 


( 16 ) 


6 








Then g can be expressed in the form (151 with m = 2 n and 


{r\, ■ ■ ■, r 2n } = {0,1 ,..., n - 1, r, r + 1,..., r + n - 1} . 
Further, we have g(x ) = x r hi(x) — h 2 (x), where 


hi(x) = Y c i n(* " Pj )’ h= E Cip i IK* “ ft)' 

*=1 o¥=i i=1 j¥=i 


Both h\ and /i 2 are Lagrange interpolation polynomials of degree at most n — 1. 
Since not all Cj are zero, neither of these polynomials is identically zero. So, if 
r ^ 1 , 2 ,..., n — 1 , then g is not the zero function. 

Hence the function g defined by (16) has at most 2n — 1 zeros in (0,oo). Of 
these, n zeros occur at x = pj, 1 < j < n. So / has at most n— 1 zeros in (0, 00 ). 
□ 


Corollary 2.3. Let r be a positive real number different from 1, 2,..., n — 1. 
Then the matrix L r defined in is nonsingular. 


Proof The matrix L r is singular if and only if there exists a nonzero vector 
c = (ci,..., c n ) such that L r (c) = 0. In other words there exist real numbers 
Ci...., c„, not all zero, such that 


n 


E 


Pi - Pi 


Pi - Pj 


= 0 


for i = 1,2,...,n. But then the function f(x) in (14) would have n zeros, viz., 
x = pi,... ,p n . That is not possible. □ 


We have proved Part (i) of Theorem o Part (iv) follows from this. If the 
inertia of L r were to change at some point r 0 > n — 1, then one of the eigenvalues 
has to change sign at r q. This is ruled out as L r is nonsingular for all r > n — 1. 


Our argument shows that if pi < p 2 < ■ ■ ■ < p n and q± < q 2 


two n-tuples of positive real numbers, then the matrix 


Pi - 9 , 


Pi-Qj 


< ■■■ < q n are 
is nonsingular 


for every positive r different from 1,2 ,,n — 1 . 

An n x n real matrix A is said to be strictly sign-regular (SSR for short) if 
for every 1 < k < n, all k x k sub-determinants of A are nonzero and have the 
same sign. If this is true for every 1 < k < r for some r < n, then we say that A 
is in the class SSR r . Sign-regular matrices and kernels are studied extensively 
in [H]. 


We have noted above that if r is any positive real number and k is any 


Pi-Qj 


Pi-Qj 

n — 1 . 


positive integer not greater than r, then every k x k matrix of the form 
is nonsingular. Let L r be an n x n Loewner matrix. Let r 7 ^ 1,2,.. 

Using a homotopy argument one can see that all k x k sub-determinants of L r are 
nonzero and have the same sign. Thus L r is an SSR matrix. If r = 1, 2,..., n— 1, 
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then the same argument shows that for k < r all k x k sub-determinants of L r 
are nonzero and have the same sign. In other words, L r is an SSRr matrix. 

Let A be any matrix with eigenvalues Ai, A 2 , • • • , A„ arranged so that |Ai| > 
| A 2 I > ■ ■ • > |A„|. The Perron theorem tells us that if A is entry wise positive, 
then Ai > 0 and Ai is a simple eigenvalue of A. (See[I3|, p. 526). Applying this 
to successive exterior powers, we see that all eigenvalues of an SSR matrix are 
simple, and the r nonzero eigenvalues of an SSRr matrix of rank r are simple. 
This proves Part (v) of Theorem o 

We now turn to proving Part (iii). Using the identity 

P[ -P r j _ Pr\p*-Pj)+PM~ 2 -P r f 2 )Pj + (Pi-PiWf 1 
Pi - Pj Pi - Pj 

we see that for every r £ R, 

L r = D r ~ 1 E + DL r _ 2 D + ED r ~ 1 , (17) 


where D is the diagonal matrix in Q and E is the n x n matrix with all its 
entries equal to one. 

By Loewner’s Theorem L r is positive definite for 0 < r < 1, and because of 
(10) it is negative definite for —1 < r < 0. Now suppose 1 < r < 2. Let x be any 
nonzero vector in the space "Hi defined in ([6). Note that this (n— l)-dimensional 
space is the kernel of the matrix E. Using o we have 

(x, L r x ) = (x, D r ~ 1 Ex ) + (x, DL r _ 2 Dx ) + (x, ED r ~ 1 x). 


The first and the third term on the right hand side are zero because Ex = 0. 
So, 

(x, L r x) = (y,L r _ 2 y), 

where y = Dx. The last inner product is negative because L r _ 2 < 0. Thus 
(x, L r x) < 0 for all x £ 'R\. In other words, L r is conditionally negative definite 
if 1 < r < 2. The same argument shows that L r is conditionally positive definite 
if 2 < r < 3 (because in this case L r _ 2 is positive definite). This was proved in 
[7] by more elaborate arguments. In particular, we have 

In L r = (l,0,n— 1), if 1 < r < 2, (18) 

and 

In L r = (n — 1,0,1), if 2 <r< 3. (19) 

We note here that if n = 3, then because of Part (iv) already proved we have 
In L r = (2, 0,1) for all r > 2. So the theorem is completely proved for n = 3. 
Let n > 3 and suppose 3 < r < 4. Now consider the space 

% = jx : y> = 0, y^pjXj = 0} 

= {x : Ex = 0, EDx = 0} . 
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This space is of dimension n — 2, being the orthogonal complement of the span 
of the vectors e = (1,1,... 1) and p = {p\,p 2 , ■ ■ ■ ,Pn)- Let x £ H 2 - Again using 
the relation we see that 

(x,L r x) = (y,L r - 2 y), 

where y = Dx. Since EDx = 0. y is in Hi, and since 1 < r — 2 < 2, we have 
(x, L r x) < 0. This is true for all x £ H 2 . So, by the minmax principle L r has at 
least n — 2 negative eigenvalues. The case n = 3 of the theorem already proved 
shows that L r has a 3 x 3 principal submatrix with two positive eigenvalues. So, 
by Cauchy’s interlacing principle, L r has at least two positive eigenvalues. Thus 
L r has exactly two positive and n — 2 negative eigenvalues. In other words, 

In L r = (2, 0, n - 2) for 3 < r < 4. (20) 

At this stage note that the Theorem is completely proved for n = 4. Now 
let n > 4, and consider the case 4 < r < 5. Arguing as before ( x,L r x ) > 0 for 
all x £ 0.2- So L r has at least n — 2 positive eigenvalues. It also has a 4 x 4 
principal submatrix with two negative eigenvalues. Hence 

In L r = (n- 2,0,2) for 4 < r < 5. (21) 

The argument can be continued, introducing the space 

h 3 = jz: = o,y^PiXj = o ,y^pfxj = oj 

= {x : Ex = 0, EDx = 0, ED 2 x = 0} 

at the next stage. Using this we can prove statement (iii) for 5 < r < 6 and 
6 < r < 7. It is clear now how to complete the proof. 

All parts of Theorem o have now been established. □ 


We end this section with a few questions. 


1. Let f(z) be the complex function defined as 


f{z) = det 


' Pi ~ Pj ' 
.Pi ~Pj. 


Our analysis has shown that / has zeros at z = 0, ±1, ±2,..., ±n — 1; 
these zeros have multiplicities n,n — 1 ,..., 1, respectively; and these are 
the only real zeros of /. It might be of interest to find what other zeros / 
has in the complex plane. 


2. When n = 3, calculations show that 


det L 3 = -(pi -p 2 ) 2 (pi -P 3 ) 2 {P 2 ~p 3 ) 2 , 
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and 


-2(pi -p 2 ) 2 (Pi ~P3) 2 (P2 -P3 ) 2 

{(pi + P2 + P 3 ){PlP 2 + PlP 3 + P2P3) + P1P2P3} ■ 


det L 4 = 


It might be of interest to find formulas for the determinants of the matrices 
L m for integers m. 

3. Two of the authors have studied the matrix P r = [(pi +Pj) r ] in [5]. It 
turns out that In P r = In L r+ i for all r > 0. Why should this be so, and 
are there other interesting connections between these two matrix families? 
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